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ABSTRACT 

Gravitational waves generated by the final merger of double neutron star (DNS) binary 
systems are a key target for the gravitational wave (GW) interferometric detectors, 
such as LIGO, and the next generation detectors, Advanced LIGO. The cumulative 
GW signal from DNS mergers in interferometric data will manifest as 'geometrical 
noise': a non-continuous stochastic background with a unique statistical signature 
dominated by the spatial and temporal distribution of the sources. Because geometri- 
cal noise is highly non-Gaussian, it could potentially be used to identify the presence 
of a stochastic GW background from DNS mergers. We demonstrate this by fitting 
to a simulated distribution of transients using a model for the DNS merger rate and 
idealized Gaussian detector noise. Using the cosmological 'probability event horizon' 
concept and recent bounds for the Galactic DNS merger rate, we calculate the evolu- 
tion of the detectability of DNS mergers with observation time. For Advanced LIGO 
sensitivities and a detection threshold assuming optimal filtering, there is a 95% proba- 
bility that a minimum of one DNS merger signal will be detectable from the ensemble 
of events comprising the stochastic background during 12-211 days of observation. 
For initial LIGO sensitivities, we identify an interesting regime where there is a 95% 
probability that at least one DNS merger with signal-to-noise ratio > unity will occur 
during 4-68 days of observation. We propose that there exists an intermediate detec- 
tion regime with pre-filtered signal-noise-ratio less than unity, where the DNS merger 
rate is high enough that the geometrical signature could be identified in interferometer 
data. 

Key words: gravitational waves - binaries: neutron stars - pulsars: general - super- 
novae. 



1 INTRODUCTION 

Three long-baseline laser interferometer GW detectors have 
been, or are nearly, constructed. The US LIGO (Laser Inter- 
ferometer Gravitational-wave Observatory) has started ob- 
servation with two 4-km arm detectors situated at Hanford, 
Washington, and Livingston, Louisiana; the Hanford detec- 
tor also contains a 2-km interferometer. The Italian/French 
VIRGO project is commissioning a 3-km baseline instrument 
at Cascina, near Pisa. There are detectors being developed 
at Hannover (the German/British GEO project with a 600- 
m baseline, which had its first test runs in 2002) and near 
Perth (the Australian International Gravitational Observa- 
tory, AIGO, initially with an 80-m baseline). A detector at 
Tokyo (TAMA, 300-m baseline) has been in operation since 
2001. The astrophysical detection rates are expected to be 
low for the current interferometers, such as 'Initial LIGO', 
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but second-generation observatories with high optical power 
are in the early stages of development; these Advanced' in- 
terferometers have target sensitivities that are predicted to 
provide a practical detection rate. 

Double neutron star (DNS) binary mergers are poten- 
tially among the strongest GW sources an d are a key search 
targe t for interferometric observatories dThorne fc Cutlet 
2002). Eight close DNS binary systems are known to ex- 
ist in the Galaxy as a subset of the observed radio bi- 
nary pulsar population ( H ulse fc Taylor] |l975|: Wolszczan 
^M). Energy loss from GW emission IITavlor fc Weisberd 
1989) causes an orbital in-spiral until the binary system 
coalesces, resulting in a burst of GWs usually described 
as a 'chirp' signal. The discovery of the DNS binary sys- 
tems containi ng PSRs J0737-3039 feurgav et al]l2003ft and 
J1756-2251 iFaulkner et alJl2005ll using the Parkes radio 
telescope brings the number of known DNS systems in the 
Galactic disk to merge within the Hubble time to four; this 
excludes the PSR B2127+11C system as it is in a globular 
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cluster and probably not formed by binary evolution (see 
Table 1). With an orbital period of only 2.45 hr, the dou- 
ble pulsar system J0737-3039A,B will coalesce in only 87 
Myr, a factor of 3.5 shorter than the coalescence time of 
the Hulse-Taylor PSR B1913+16 system, which yielded the 
first evidence for gravitational radiation. Predicted merger 
rates are dominated by the J0737-3039 system, because of 
its proximi ty, short coalescence time and the difficulty of its 
discovery (|Faulkner et al.ll2005l) . In contrast, PSR J1756- 
2251's parameters are similar to those of previously known 
pulsars; thus, its very recent addition isn't expected to sig- 
nificantly alter the predicted merger rate. 

In order to estimate the DNS coalescence rate, Kalogera 
et al. (2001) used a semi-empirical approach, based on the 
physical properties of known DNS systems and pulsar survey 
selection effects, to obtain scale factors that correct for the 
unobserved fraction of existing systems. In an alternative 
approach, Kim et al. (2004) generated the probability dis- 
tribution of merger rates using Monte Carlo simulations in 
which they populate a model galaxy with a specified number 
of NS binaries whose properties are those of known systems; 
they also model survey selection effects to deduce detection 
rates and their statistical significance. With this model and 
more up-to-date pulsar sur vey results that include J0737- 
3039, IKaloeera et alJ (120041) presented updated bounds for 
the coalescence rate for Galactic disk DNS systems, yield- 

and derived detection rates for 
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jPiranl \im& lAndol l2Q04t) assume a simple 
parametrization of merger rates with merger probability 
P m (t) cx t a , where t is the time from formation of the 
binary and a is a constant; the distribution of mergers 
as a function of time is simply given by the convolution 
of the star formation rate with the distribution P m (t). 
iReeimbau et alJ (120051) . using numerical simulations, derive 
a probability distribution for DNS merger time, r, of the 
form P(t) — 0.087/r, and find that mergers are possible 
between 2 x 10 yr and the age of t he Universe. Using both 
evolutionary and statistical models. IReeimbau et alJ (I2005T) 
find a Galactic merger rate of 1.7 x 10~ 5 yr -1 , sim ilar to the 
lower bound calculated bv IKaloeera et all ll2004h . 

A purely theoretical approach to estimate coalescence 
rates can be based on the predictions of population synthesis 
models that attempt t o simulate the entire pu lsar population 
jLipunov et alJll997t iBelczvnski et al.l F2002i. However, the 
size of the parameter space and associated uncertainties in- 
volved seriously limit this way of inferring merger rates. For 
instance, new theoretical models of merger formation can 
have considerable impact: Recent simulations by Chaurasia 
and Bailes (2005) show that significant kick velocities im- 
parted to neutron stars at birth give rise to highly eccentric 
orbits resulting in accelerated orbital decay and greatly re- 
duced lifetimes. Their work indicates that a non-negligible 
fraction of DNS systems could merge within times as short 
as 10 yrs. This implies that there could be a selection ef- 
fect favouring the observation of those binary pulsar systems 
that coalesce over relatively long times. 

Our interest is in the detection of a GW background 
generated by DNS merger events throughout the cosmos. 
As a starting point, we adopt the Galactic coalescence rates 
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Figure 1. The differential DNS merger rate as a function of red- 
shift 2 using the lower Galactic rate of 17 Myr -1 from Kalogera 
et al. (2004) and four merger times, r = (0,1,2,5) Gyr, assumed 
the same for all systems. The curves are calculated using the SFR 
evolution factor obtained from the parametrized model labelled 
SF2 in Porciani & Madau (2001) and a flat-A cosmology with 
f^m = 0.3 and Q/± = 0.7, assuming a SFR cut-off at z = 5; the 
cut-off projects forward in time to the end points of the time- 
delayed curves. A time delay of 1—5 Gyr has minimal effect on 
di?DNS /da at small z compared to the differential rate locked to 
the evolving SFR. 



of IKaloeera et all (|2C)04) as a basis for estimating the rate 
density in our local intergalactic neighbourhood. To project 
this throughout the Universe, we assume that the DNS sys- 
tem formation rate follows the evolving star formation rate 
(SFR), for DNS systems with relatively short merger times 
of tens to hundreds of millions of years. In this way, we em- 
ploy the updated Galactic DNS merger rate estimates as in- 
put to a model that describes the observer's evolving record 
of DNS mergers in terms of a 'probability event horizon' 
(PEH). We describe how the detectability of DNS mergers 

evolves with observation time. 

We use the PEH concept dCoward fc Burmanll2005l) to 
describe the cumulative GW signal from cosmic DNS merg- 
ers and show that the signal is dominated by the spatial and 
temporal distribution of events. Instead of describing such 
a signal as 'stochastic', we argue that the signal, although 
temporally random, has a unique signature. We use the term 
'geometrical noise' to distinguish it from the stationary noise 
usually associated with a primordial GW stochastic back- 
ground. 



2 THE COSMOLOGICAL PROBABILITY 
EVENT HORIZON 

In the standard Friedmann cosmologies, one can express the 
differential DNS merger rate as the merger rate in the red- 
shift shell 2 to z + dz: 



di?DNS 



dV roe(z) 
dz 1 + z 



dz , 



(1) 



where dV is the cosmology-dependent co-moving volume ele- 
ment and -Rdns(z) is the all-sky (4-7T solid angle) DNS merger 
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Table 1. DNS systems containing radio pulsars arranged in order of increasing estimated coalescence 
time. The first entry is the double pulsar system; the second system is in a globular cluster (M 15); 
the third is the Hulse- Taylor binary. Here, e is the orbital eccentricity, t c is the pulsar characteristic 
age and tqw i s the time remaining to coalescence due to emission of gravitational radiation. The 
relationship between a pulsar's characteristic and true ages is not clear, but we suppose the total 
coalescence time from birth to be t c + tq^i . Bold type indicates DNS systems that will coalesce in 
less than 10 10 yr; the gobular-cluster system B2127+11C is excluded because its formation history is 
probably different from those of the Galactic-disk systems: it is likely that this system has undergone 
partner exchange from close interactions. References: (1) Lyne et al. (2004); (2) Burgay et al. (2003); 
(3) Anderson et al. (1990); (4) Hulse & Taylor (1975); (5) Faulkner et al. (2005); (6) Wolszczan 
(1990); (7) Lyne et al. (2000); (8) Nice et al. (1996); (9) Champion et al. (2004). Adapted from 
Faulkner et al. (2005) 
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(ms) (hr) (M ) (Myr) (Myr) 



J0737 3039A,B 


23, 2773 


2.45 


0.088 


2.58 


210,50 


87 


(1,2) 


B2127+11C 


31 


8.04 


0.681 


2.71 


969 


220 


(3) 


B1913+16 


59 


7.75 


0.617 


2.83 


108 


310 


(4) 


J1756-2251 


28 


7.67 


0.181 


2.57 


444 


1690 


(5) 


B1534+12 


38 


10.10 


0.274 


2.75 


248 


2690 


((>) 


J1811-1736 


104 


449 


0.828 


2.6 


900 


1.7 x 10 6 


(7) 


J1518+4904 


41 
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0.249 


2.61 


16,000 


2.4 x 10 6 
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J1829+2456 


41 


28.2 


0.139 


2.5 




6 x 10 4 
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rate, as observed in our local frame, for sources out to red- 
shift z. Source rate density evolution is accounted for by 
the dimensionless evolution factor e(z); this is normalized 
to unity in the present-epoch universe (z = 0) so ro is the 
z — rate density. The Galactic DNS merger rate is con- 
verted to a rate per uni t volu me, ro, using the conversion 
factor 1CP 2 from lAndol l)2004l ) for the number density of 
galaxies in units of Mpc -3 . The (1 + z) factor accounts for 
the time dilation of the observed rate by cosmic expansion, 
converting a source-count equation to an event rate equa- 
tion. 

We assume a 'flat-A' cosmology, in which the density 
parameters of matter and vacuum energy sum to unity; 
we use O m = 0.3 and Q,a = 0.7 for their present-epoch 
values and take Ho = 70 km s _1 Mpc -1 for the Hubble 
parameter at z = 0. With the assumption that the DNS 
system formation rate follows the evolving star formation 
rate, e(z) is obtained by normalizing a SFR evolution model 
to its z = value. Presently there is no consensus among 
astronomers on a model of star formation history at high 
z, so we shall use the pa rametrized model labelled SF2 of 
Pgrcigii^^Mg^3 {200 1 ) scaled to the flat-A cosmology; see 
Coward fc Burmad pOO^) for details. The cumulative DNS 
merger rate .Rdns(z) is calculated by integrating the differ- 
ential rate from the present epoch to redshift z: 

Rbns(z)= f (di? DNS /dz) dz. (2) 
Jo 

Whereas core-collapse supernova events (because of 
their short-lived progenitors) closely follow the evolving star 
formation rate, DNS merger events are generally signifi- 
cantly delayed with respect to star formation. In order to 
test the sensitivity of dffoNs/dz to this delay, we have cal- 
culated it under the simplified assumption that all such 
events are delayed by the same merger time r. Figure Q 
plots d7?DNs/dz for four choices of r, namely 0, 1, 2 and 5 



Gyr. The three delayed curves were obtained by converting 
e(z) to a function of cosmic time (or, equivalently, look-back 
time), shifting it by r and then converting back to a function 
of z. Cosmic time and z are related by the evolving Hub- 
ble parameter H(z), which can be expressed in terms of th e 
contributions of matter and vacuum energy dPeebleslll993Tl : 

h{z) = H(z)/H = [« m (l + zf + Q A ] 1/2 (3) 

for a flat-A cosmology, and cosmic time is expressed by the 
formula: 

T cos (z)= f [(l + z)ft(z)]-Mz. (4) 
Jo 

Figure 1 shows that a 1-2 Gyr merger time does not 
significantly alter di?DNs/dz, especially for small z. Increas- 
ing t has the effect of reducing the available volume for 
observing DNS mergers and, for r > 5 Gyr, coalescences are 
restricted to the volum e bounded by z sa 1. S o, for systems 
such as those used by iKaloeera et alJ i2004j) in their rate 
estimates, we can safely assume that they quite closely fol- 
low the evolving star formation rate. For definiteness, in the 
following calculations, we assume that the DNS merger rate 
approximately follows the SFR evolution factor e(z). 

As the DNS mergers throughout the Universe are inde- 
pendent of each other, their distribution is a Poisson process 
in time: the probability e for at least one event to occur in 
a volume out to redshift z during observation time T at a 
mean rate Runs(z) is given by an exponential distribution: 

p(rO l;7? DNS (z),T) = l-e- flDNs(z)T = e, (5) 

with mean number of events jV £ = Runs( z )T — |ln(l — e)|. 
The PEH is defined by the minimum distance, or redshift 
zf EH , for at least one event to occur over some observation 
time T, with probability above some selected threshold e. We 
find zf EH (T) by fixing e and solving the above condition nu- 
merically, thus defining the PEH. In practice, we shall take 
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Table 2. The probability for at least one DNS merger to occur 
within the volume bounded by the luminosity distance = 20 
Mpc for one year of observation (Initial LIGO) and by d^ = 200 
Mpc for one day (Advanced LIGO) assuming the Galactic DNS 
merger rate limits, -RdnSi from Kalogera et al. (2004). 



Galactic -Rdns 


p(d L = 20 Mpc; 1 yr) 


p(d L = 200 Mpc; 1 day) 


(Myr" 1 ) 


Initial LIGO 


Advanced LIGO 


292 


0.13 


0.86 


17 


0.008 


0.11 



p= 11.7 



1*1 ( ruao \ 
2.56 V r J 



(6) 



where tligo = 20 and 200 Mpc for Initial and Advanced 
LIGO respectively and l*^ 2 , the angular sensitivity function 
of the detectors, averages to 2.56 over all sky and source 
positions. By substituting d^f H (T) for r in equation ©, 
we obtain the DNS merger probability event horizon of the 
optimal signal-to-noise ratio: 



Pe (T) 



11.7 



?"LIGO 



2.56 \dlf H (T)J 



(7) 



Assuming p ~ 11.7, corresponding to the detection 
threshold for a DNS mer ger in Advanc ed LIGO for an 'av- 
erage' source at 200 Mpc l|Hughesll2002fl . one can model the 
evolution of p as a function of duty cycle or observation 
time. Figure |3] plots Po™(T), showing the PEH for at least 
one DNS merger to occur during observation time T, with 
P Po^k(T), for the upper and lower uncertainty limits of 
the local DNS rate densities iKaloeera et al 1l2004l) . For Ad- 
vanced LIGO, at least one DNS merger with a p > 11.7 will 
be detectable using a matched filter every 12-211 days with 



10 



B 1 



e = 0.95, corresponding to a 95% probability of observing at 
least one event within 2, so N c — 3 is the mean number of 
events. Converting zf EH (T) to luminosity distance d?[f H (T), 
and differentiating with respect to T then yields vf BH (T): 
the PEH velocity, which describes the rate at which obser- 
vations penetrate into the low-probability 'tail' of the event 
distribution. 

Table [5] uses equation (J^J to calculate the probabil- 
ity for at least one DNS merger to occur within the vol- 
ume bounded by — 20 Mpc for one year of observa- 
tion (Initial LIGO) and within d^ = 200 Mpc for one day 
(Advanced LIGO) using the maximum and minimum val- 
ues of_tLe_Cj^la£tic^DNS merger rate Rons as estimated 
bv iKaloeera et al.l <2004) . namely 292 and 17 Myr" 1 . The 
two distances, 20 and 200 Mpc, correspond approximately to 
the detection thresholds for DNS merger detection for Initial 
and Advanced LIGO respectively. Figure |5] plots d[ £ EH (T), 
using the upper and lower Galactic merger rate limits from 
Kalogera et al. (2004). 



3 THE PEH FOR LIGO TO DETECT DNS 
MERGERS 

For a DNS merger at distance r from the detector, the opti- 
m al GW signal-t o-noise ratio, p, assuming matched filtering 
is (iHuehesfeool) : 
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Figure 2. The PEH evolution for DNS mergers as a function of 
observation time assuming the Galactic merger rate limits from 
Kalogera et al. (2004): 292 Myr -1 (dashed curve) and 17 Myr -1 
(solid line) , and merger rate evolution following the instantaneous 
SFR. The horizontal lines show the 200 and 20 Mpc horizons, 
indicating the sensitivity limits for DNS coalescence detection 
assuming optimal filtering for Advanced and Initial LIGO. 
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Figure 3. The 95% PEH for the temporal evolution of the signal- 
to- noise ratio p^ EH (T) for Advanced and Initial LIGO using the 
local Galactic rate limits for DNS mergers from Kalogera et al. 
(2004): 292 Myr" 1 (dashed curve) and 17 Myr" 1 (solid line). For 
Advanced LIGO with detection threshold p > 11.7 (assuming 
optimal filtering), at least one DNS merger will be detectable 
during 12—211 days of observation. For initial LIGO sensitivities, 
DNS mergers will be detectable in 22—380 yrs of observation and, 
for a p ~ 1, in 4—68 days. In the region below p = 11.7, individual 
mergers are unlikely to be detected, but the cumulative signal 
will manifest as 'geometrical noise' in interferometric data. 



a 95% probability. For initial LIGO sensitivities, events with 
the same p will be detectable at a rate of one per 22-380 yr 
and for p ~ 1, one every 4-68 days. Although DNS mergers 
with p < 11.7 will be present in the data, they are unlikely 
to be temporally resolved. 

Here, we consider a signal comprised of an ensemble of 
chirp signals that result from the final mergers of cosmo- 
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logical DNSs. Even though the individual signals will not 
overlap, the cumulative background is stochastic because 
one cannot predict with certainty the exact temporal evo- 
lution of the signal. But, the GW amplitude distribution 
will be very different from a continuous astrophysical GW 
background. For a signal that consists of a sum of many 
continuously emitting, randomly distributed sources, such 
as the GW signal from the Galactic population of white 
dwarf-white dwarf or white dwarf-NS binaries, the central 
limit theorem predicts that the GW amplitude distribution 
will be Gaussian. The term 'confusion noise' is often used 
to describe the signal from the Galactic population of white 
dwarf- white dwarf binaries. In contrast, the GW background 
from DNS mergers has a non-Gaussian signature that is 
characterized by the spatial distribution of the sources, so 
we use the term geometrical noise to describe it. 



THE PEH TECHNIQUE APPLIED TO 
SIMULATED DATA 



T3 



ICoward. Burman fc Blair] <l2002a.bl) developed a procedure 
to simulate time series of the GW signals from supernovae 
occurring throughout the Universe. The method has pro- 
vided a tool to probe the statistical signature of cosmological 
'standard candle' GW sources that follow the evolving star 
formation rate. In its simplest form, the simulation outputs a 
distribution of events in redshift as a function of observation 
time. Here, we use this procedure to extract the PEH sig- 
nature from a distribution of observer- source distances and 
observation times. The PEH technique is a procedure that 
records the time and distance of successively closer events, 
(ti, n), that satisfy the condition r^+i < re each event in the 
PEH time series occurs at a shorter distance than the pre- 
ceding one. A key feature of the procedure is that it probes 
the rate of an observer's penetration into the small proba- 
bilit y 'tail' of an event distribu tion. 

ICoward fc Burmanl j2005l) applied this technique to the 
observed gamma-ray burst (GRB) redshift and observation- 
time distribution and found that it can be used to set limits 
on the local rate density of GRBs. Because it is sensitive 
to the distribution of events in the small probability tail, in 
this case corresponding to smaller volumes, the procedure 
'tracks' the temporal evolution of close events as detected 
by an observer. If a time series is comprised of more than 
one class of event, then the PEH technique can identify these 
distributions b ased on the distribution o f PEH data. This 
was shown by ICoward fc Bur man ( 2005J) , where the PEH 
model fitted to the data indicated that the anomalously close 
GRB 980245 probably does not belong to the 'classical' GRB 
population. 

As a further proof of concept, we apply the PEH tech- 
nique to a distribution of GW amplitudes from a simu- 
lated distribution of DNS mergers throughout the Universe. 
We use a DNS dimensio nless characte ristic GW amplitude 
h c = 4x 1CT 21 at 20 Mpc <Schutzll99 lb and a non-redshifted 
characteristic frequency f c = 200 Hz. The results are plotted 
in figure 2] The 5% and 95% PEH curves are shown along 
with the output from the PEH simulation shown as data 
points. The 5% PEH curve is interpreted as the 'null PEH', 
as there is a 95% probability that no events will be observed 
at an amplitude above this threshold. 
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Figure 4. The PEH technique applied to a distribution of GW 
amplitudes from simulated DNS mergers throughout the Uni- 
verse assuming the upper Galactic rate limit of 292 Myr - 1 from 
Kalogera et al. (2004). We use a DNS characteristic GW ampli- 
tude he = 4x 10 -21 at 20 Mpc (Schutz 1991) and a non-redshifted 
characteristic frequency f c = 200 Hz. The 5% (or 95% probability 
of no events occurring) and 95% PEH curves for DNS mergers are 
shown, along with the output of the PEH simulation applied to 
a distribution of characteristic GW amplitudes modelled on DNS 
mergers, shown as filled circles. Also shown are the 5% and 95% 
'Gaussian- noise PEHs', modelled on idealized interferometric de- 
tector noise-see text for details. The figure shows that the DNS 
merger PEHs evolve much faster than the Gaussian-noise PEHs. 



In order to develop the PEH concept into a technique 
for GW data analysis, we need to incorporate a method for 
treating detector noise. For 'standard candle' sources, this 
can be done by working in terms of progressively increasing 
amplitudes rather than progressively decreasing redshifts, 
as shown in figure 5. Noise transients of varying amplitudes 
mimic sources at varying redshifts, enabling the introduc- 
tion of a 'noise PEH'. Here we use a 'Gaussian-noise PEH'; 
this is based on a standard zero-mean Gaussia n distribution 
with standard deviation s d — h rms ^J f s /2f c , modelled on 
idealized interferometric detector noise with h rma denoting 
the root-mean-squa re amplitude of the n oise and f s the sam- 
pling freauencv-see lArnaud et alJ (|200.?f l . Figure 5 shows the 
5% and 95% Gaussian-noise PEHs. 

It is clear from the small gradient of the Gaussian-noise 
PEH curves that the temporal evolution of transients, origi- 
nating from the tail of the distribution is slow. The striking 
feature of this figure is the gradient of the PEH curves for 
the source amplitude distribution: it is much bigger than for 
the Gaussian-noise case, because of the distance-amplitude 
relation. We define this type of PEH data as 'geometrical', 
because it is linked to the spatial distribution of sources. 

We note that although we have incorporated SFR evolu- 
tion and a cosmological model to calculate the DNS merger 
rate at high z, the PEH will be most sensitive to events that 
are comparable to or greater than the detection threshold 
of the detector. For Advanced LIGO, this will be of order 
several hundred Mpc (z » 0.05), implying that the PEH 
will be dominated by sources in the local Universe where 
cosmological effects are minor. 

Given that a geometrical signature is present in a time 
series, one can fit a PEH model to the data, with the local 
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rate density as a free parameter. In the presence of noise, this 
signature will be modified, but if the detector noise is well 
characterized, a fitting function constructed from the sum 
of both noise and geometrical distributions can be used. For 
small signal-noise ratios, the noise will dominate the PEH 
data for short observation times. But as observation time in- 
creases, the geometrical signature will grow at a faster rate 
than the noise. It is possible that the GW signal from DNS 
mergers can be identified before that of a single merger, 
occurring above the detector threshold, as shown in figure 
[3] We propose that there exists an intermediate detection 
regime, with pre-filtered signal-noise-ratio of less than unity, 
where the DNS merger rate is high enough that the geomet- 
rical signature could be identified in interferometer data. 
This type of search provides an opportunity to set upper 
limits on the DNS merger rate. 



5 DISCUSSION 

The cumulative signal from transient GW sources at cos- 
mological distances is commonly described as a stochastic 
background because of the temporal randomness of the in- 
dividual events. But we have shown here that, for a cos- 
mic population of DNS mergers, the GW cumulative signal 
will be dominated by the spatial distribution and tempo- 
ral evolution of the sources. We use the term 'geometrical' 
noise, because the GW amplitude distribution composed of 
DNS mergers is dominated by the spatial distribution of the 
sources. 

Interestingly, geometrical noise from cosmic DNS merg- 
ers (and similar transient GW sources) is linked to the sensi- 
tivity of the detector and the cumulative event rate. If a GW 
detector could probe DNS mergers to distances correspond- 
ing to the time when they first occurred in the Universe, the 
detector could potentially resolve all events and there would 
be no geometrical noise from unresolved signals. However, 
if the cumulative rate measured in our frame were to be so 
high that the individual chirps would overlap, correspond- 
ing to a signal of duty cycle greater than unity, then the 
resulting signals would manifest as geometrical noise. This 
is unlikely for DNS mergers: the individual events are not 
expected to overlap but to be nonetheless not individually 
detectable because of the small signal-to-noise ratio. 

We highlight the concept that the detectability of a GW 
stochastic background from DNS mergers throughout the 
Universe and of individual DNS mergers is linked to the ob- 
server via the instrument duty cycle, the cumulative event 
rate and the sensitivity of the detector. The PEH model is 
applied above to the detectability of DNS mergers with Ad- 
vanced and Initial LIGO; it shows the temporal evolution of 
the GW signal from many faint unresolved events converg- 
ing to detectable single events of high signal-to-noise ratio. 
Detection strategies based on the PEH model could utilize 
the geometrical signature of the temporal evolution of the 
chirp GW amplitude distribution. We plan to investigate the 
PEH technique further using simulated GW background and 
non- Gaussian noise data. 

The PEH method could also be used to set constraints 
on the black hole-black hole (BH-BH) merger rate, which 
is presently highly uncertain. Although occurring at a lower 
intrinsic rate than DNS mergers, the GW emission from BH- 



BH mergers is expected to be significantly higher, imply- 
ing an enhanced detection rate because of the larger detec- 
tion volume. Estimates range from 30-4000 yr _1 assuming 
a BH mass of 10 Mr and a detection horizon of z — 0.4 
jThorne fc Cutlejl2002|) for Advanced LIGO sensitivity. We 
plan to use the BH-BH mass and merger rate distributions 
to determine if a PEH model could be developed and ap- 
plied to interferometer data to set limits on the rate and 
mass distributions. 

Although we have shown how the cumulative GW signal 
from DNS mergers should manifest in a time series, we have 
not yet tested how robust the above procedure is to the 
non-stationarity of the data. For instance, if a time series is 
plagued with non-Gaussian transients, then the effectiveness 
of the fitting procedure will be reduced, unless a good model 
for the transients is available. 

We assume optimal (matched) filtering can be applied 
to the candidate transients, but the filters that will be most 
effective will not be exactly matched to individual events. 
The 'norm' and 'mean' filters, although not optimal, should 
perform at a significant fraction of a matched filter for 
transients where th e signal is not completely characterized 
iArnaud et all2003Tl . We emphasize that using the PEH in a 
signal processing context is based on fitting to outliers using 
a model distribution and not 'detecting' individual events. 
Nonetheless, it is important to filter the data to select as 
many transients that are candidates for DNS mergers as 
possible. As a test of the robustness of the procedure using 
different filters, we plan to inject simulated cosmologically 
distributed GW transients into non-Gaussian noise. 

We note that the simulated GW amplitude distribution 
will be modified by the detector antenna pattern, which we 
have not included in our simulations, where we used the 
mean value. In applying a PEH model to real data, account 
must be taken of the modulation and distortion of the signal 
by detector characteristics and by both detector and envi- 
ronmental noise. We plan to extend the current simulations 
by incorporating such effects as part of a program to test 
the sensitivity of the PEH fitting procedure. 
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